The Density of States in High-T c Superconductors Vortices 
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We calculated the electronic structure of a vortex in a pseudogapped superconductor within a 
model featuring strong correlations. With increasing strength of the correlations, the BCS core 
states are suppressed and the spectra in and outside the core become similar. If the correlations 
are short-range, we find new core states in agreement with the observations in YBa2Cu3C>7_i and 
Bi2Sr2CaCu2 08+5 • Our results point to a common phenomenology for these two systems and 
indicate that normal-state correlations survive below T c without taking part in the overall phase 
coherence. 

PACS numbers: 74.20.-z, 74.25.Jb, 74.60.Ec 
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Among the anomalous properties of high-temperature 
superconductors (HTS), one of the most intriguing is the 
pseudogap phenomenon observed in the normal state of 
underdoped systems. The pseudogap state is marked 
by a substantial decrease of the one-particle density of 
states near the Fermi energy, and can be observed in 
several spectroscopic, thermodynamic, and transport ex- 
periments (see Ref. |l| for a review). Its elucidation is 
widely thought of as one of the keys to the understanding 
of high-temperature superconductivity. Unfortunately, 
detailed studies of the low-lying excitations in the pseu- 
dogap state are impossible due to the emergence of the 
superconducting state at low temperature. Vortex-core 
spectroscopy is perhaps one way into this problem: it 
has indeed been stated, on the basis of experimental 
findings that the local density of states (LDOS) in 
vortex cores below T c reflects the properties of the (pseu- 
dogapped) normal state above T c . More precisely one can 
argue that the physical state inside vortices should reflect 
the properties of the normal state at low temperature in 
the boundary conditions set by a gradual vanishing of the 
superconducting order in the core. A detailed analysis of 
this density of states may therefore lead to useful insights 
into the nature of the pseudogap state itself. 

Some characteristics of the vortex cores seem to be 
now well established in YBa2Cu307_,5 (YBCO) and 
Bi 2 Sr 2 CaCu208+(5 (BSCCO), the two systems for which 
experimental results are available^, ||, ||, [|: the LDOS in 
the core does not display a central peak as predicted by 
the BCS theory, but particle-hole symmetric low-energy 
states; noticeably, the amplitude of the core states is 
much larger in YBCO than in BSCCO. In BSCCO, more- 
over, the core LDOS shows remarkable similarity (when 
properly thermally broadened) to the density of states 
observed above T c in the pseudogap regime. 

Theoretically, a number of calculations have been done 
in the past years and competing models of HTS lead to 
increasingly detailed predictions. The first calculations 
were undertaken within the Bogoliubov-de Gennes for- 
mulation of the BCS theory §, ff |, g For a d x 2_ y 2 
vortex the LDOS shows a broad zero-bias conductance 



peak (ZBCP) and a four- fold symmetric shape. This pic- 
ture is not consistent with the experimental observations 
in YBCO and BSCCO. In another BCS- type scenario, a 
secondary order parameter component of a different sym- 
metry is generated locally by the magnetic field and the 
resultant pair potential no longer has nodes: the central 
peak would thus be suppressed or split 0. The absence 
of ZBCP has also been tentatively explained by the prop- 
erties of the c-axis tunneling matrix elements |DJ]. Many 
workers in the field believe however that the absence of 
central peak cannot be explained in the BCS framework, 
and needs some kind of fundamental extension of the 
BCS theory that would take into account strong corre- 
lations. The first attempts were made within the SO (5) 
theory and predicted an insulating core H [l§, but the 
detailed predictions of this model seem not to agree with 
experiment. Some t-J model calculations have predicted 
a ZBCP split by an induced s-wave order parameter at 
low doping jl4| [l]|. Very recently, calculations based on 
the SU(2) slave-boson theory were reportedQ. The au- 
thors make detailed predictions for the LDOS in vortex 
cores, which do not seem to correspond to the experi- 
mental data available today and which should allow to 
test in a significant way the slave-boson theory of HTS. 



In this Letter we present calculations based on the 
Cooperon-propagator description of HTS. We show that 
this approach can explain several characteristics of the 
vortex-core spectra in YBCO and BSCCO. In addi- 
tion, it allows us to give a tentative explanation for 
one of the (many) puzzles of HTS, namely whether the 
gap and pseudogap are of similar (superconducting) ori- 
gin or whether they correspond to completely differ- 
ent physical processes. This point will be commented 
upon in the discussion below. First, we briefly review 
the Cooperon-propagator description of superconductiv- 
ity and the model used in the numerical calculations. 
In the Kadanoff-Martin theory of superconductivity |L7[ , 
which is entirely based on correlation functions, the elec- 
tronic self-energy takes the formp8[ 
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£(r,s,"r) = - Jdfds(T T {A(r,r,T)A\s 7 s,0)}) 

xg (s-f,-T), (i) 

where A(r,s,r) = F(r — s)ip^(r, r)ipi (s, r), V(r) is the 
pairing interaction, and Go(r,r) is the free propagator. 
In this formalism the superconducting order is charac- 
terized by the long-range properties of the pair correla- 
tion function. Eq. (|l|) is valid both above and below T c 
provided higher order correlation effects are not impor- 
tant, and we will assume that it remains valid regardless 
of the model or approximations involved in calculating 
(AA^). In particular, the BCS theory is recovered by 
setting (AA f ) = (A) (At) where (A) is the BCS order 
parameter. Following recent studies)!^] we model the 
correlation function (Cooperon propagator) through: 



<T T {AAt}) 



A p (r-f)R(£/g)A p (s-s) 

(A s (r-r)F(£/OA s (s-s) T > T c 



A s {r,f)A*(s,s) 



(2) 



The first term characterizes the strength and symmetry 
(A p ) and the range (g) of the normal-state correlations 
leading to the pseudogap below the temperature T*, with 
i = i|r + r — s — s\ and R(x) = e~ x . The second term 
accounts for the phase-coherence properties of the su- 
perconductor and the phase fluctuations responsible for 
destroying superconductivity. Above T c , A s must be un- 
derstood as the modulus of a local BCS-like order pa- 
rameter. The phase physics is contained in the function 
F(x) ~ e~ x which is the phase-phase correlation func- 
tion. F thus reduces to a function of the distance I with a 
correlation length £(T) diverging at T c in the Kosterlitz- 
Thouless fashion. Below T c , F is assumed to factorize, as 
in the BCS theory, into the product of two local phases. 
These phases are embodied in the amplitude A s which 
becomes position dependent and complex. A more elab- 
orate model should include dynamic correlations. We 
showed, however, that Eq. (||) leads to reasonable agree- 
ment with some aspects of the experimental density of 
states as a function of T in homogeneous sy stems [O. 
Eq. (|J) may be viewed as a non-local generalization of 
the Ansatz A 2 = A 2 + A 2 introduced in Ref. 

In what follows we concentrate on the vortex state at 
T = 0, then A s plays the role of a BCS gap and van- 
ishes in the vortex cores while its phase winds by 2-7T in 
response to the applied field. Eq. (j|) basically means 
that the Cooperon propagator has a short-range part in 
addition to the long-range part described by A s . Our 
central assumption is that these short-range normal-state 
correlations are not involved in the phase coherence and 
remain finite in the cores. This view is supported by re- 
cent NMR experiments showing that the pseudogap is 
not affected by a strong magnetic field poj. We also note 
that similar considerations apply to the spinon and holon 
fields in theories based on spin-charge separation plL E2|. 



Using Eqs. (||) and (||) we calculate the Green's func- 
tion for an isolated d x 2_ y 2 vortex on a two-dimensional 
lattice. We neglect Landau-level quantization, as ap- 
propriate for low magnetic fields. The modulus of A s 
is proportional to tanh(r/a), with a the lattice parame- 
ter. The pseudogap A p has d x 2_^i symmetry in agree- 
ment with photoemission data|23] and is assumed uni- 
form in space for simplicity. The self-energy is rewritten 
as = Ti u (u>) + where E u is uniform and Ej 

contains the inhomogeneous term. The uniform part of 
the Green's function is evaluated as Q u = (C?^ 1 — Eu) -1 in 
reciprocal space, using a 1024 x 1024 fc-point mesh. The 
band structure is taken from Ref. |24| with the Fermi level 
Ep located 10 meV above the van-Hove singularity. This 
yields zero-field DOS's in qualitative agreement with ex- 
periment. The Green's function is finally evaluated as 
Q = (Q^ 1 — Ei) _1 on a M x M real-space mesh with the 
vortex at the center. This step necessitates the inver- 
sion of two large matrices (M 2 x M 2 ) for every energy u>. 
We have checked that the calculations are well converged 
with M = 41 and reproduce the published results of the 
BCS theory when A p = [see Fig. |(c)]. 

Our main results are presented in Fig. [I]. Fig. |l](a) 
shows the vortex-core LDOS N(r = 0, oj) and the zero- 
field DOS Nq(lu) for a pseudogapped superconductor. In- 
stead of a ZBCP, we find two core states in the vortex 
spectrum. Our numerical analysis shows that the sub- 
gap peaks appear only for moderate values of A p /A s 
and low values of g. These peaks grow as A p /A s and/or 
g decrease, and eventually merge to form the broad BCS 
peak as A p — » 0. An intermediate example is displayed 



in Fig. U\(b). In fact, weak subgap structures exist also in 
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FIG. 1: Vortex-core and zero-field DOS for the parameters 
g = 5a and (a) A = ,/AjTAj = 40 meV, A p = 4A S ; (b) 
A = 20 meV, A p = 2A S . (c) BCS limit: A p = and A s = 
40 meV. Inset: Halved core- level separation as a function of 
A for g — 5a (filled symbols) and g — 10a (empty symbols). 
Circles are for A p — 2A S and triangles for A p = 4A S . No 
core states are found for g — 10a and A p = 4A S . 



the zero-field DOS: the vortex seems to induce a localiza- 
tion of preexisting low-energy states. We therefore expect 
those pseudogap-induced core states to be independent 
of magnetic field as observed in recent experiments [p5[ . 
On the contrary the BCS core states were shown to split 
at high vortex density JlCj] ; a field-dependent splitting is 
also expected in scenarios based on a secondary id xy or- 
der parameter jj| . 

It was found experimentally that the coherence peaks 
are suppressed in the core||, @, |], unlike in our results, 
Fig. |l|(a) and (b) . In our calculations the spectral weight 
in a small energy window (~ 4A) around E-p is locally 
conserved, in contrast to experiments. This discrepancy 
might be due to a vortex charging effect or to our as- 
sumption of static correlations in the model Eq. (||) . Al- 
though, at this stage, the model is unable to account for 
the observed behavior of the coherence peaks in the core, 
we tentatively argue below that it can well explain the 
properties of the core states in both YBCO and BSCCO. 

The energy AE of the core states is plotted in the inset 
of Fig. [| as a function of the gap scale (Aj; + A^) 2 . The 
calculated core-state energies increase with the gap width 
and the correlation length g but depend little on A p / A s . 
They are in good general agreement with the experimen- 
tal observations, i.e., 5.5 meV in a YBCO sample with a 
gap of 20 meVf§, 7 meV for a 32 meV-gap BSCCO §, 
and 14 meV for a 45 meV-gap BSCCO|@. A recent 
thorough investigation of the BSCCO system indicates 
a roughly linear relationship between AE and the en- 
ergy gap p5|] . Moreover, the fact that the observed peaks 
are sharper in YBCO than in BSCCO can be well ex- 
plained by assuming a smaller value of A p /A s in YBCO 
[see Fig. |l](b)] . This assumption is consistent with the ob- 
servation that optimally-doped YBCO, unlike optimally- 
doped BSCCO, seems to have no pseudogap phase above 

We now discuss the behavior of the LDOS in the vicin- 
ity of the core. N(r, lu) is displayed in Fig. ||(a) as a func- 
tion of oj for r along x , illustrating the weak localization 
of the core states. The angular average N(r, uj) plotted in 
Fig. |(b) at the core-state energy shows the exponential 
decay of the LDOS, as observed experimentally Q and 
in contrast to the algebraic decay of BCS core states. 
The calculated decay length depends on energy and is 
A/a = 6.7 for the state above E-p and 4.4 for the state 
below Ep. The average 5.6 compares well with the value 
of 5.7 ± 0.8 quoted in Ref. (using a = 3.83 A for the 
BSCCO lattice parameter). Fig. |^(b) also shows the 
absence of a significant four-fold anisotropy in N(r,Lu) 
around the core. It is known from semiclassical ap- 
proaches that a finite lifetime reduces the anisotropy P7J . 
The pseudogap correlations have the same effect in our 
calculations. As an illustration we plot in the inset the 
ratioa(r) = f Q dui[N(rn,ui) — N(rx,uj)]/f dujN(r,uj) 
which measures the relative anisotropy of the LDOS at 
subgap energies (n is the nodal direction). In the BCS 
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FIG. 2: LDOS near the vortex core for the model parame- 
ters as in Fig. |l|(a). (a) N(r,uj) at all sites between ( — 10,0) 
and (10, 0) (top to bottom). The curves are shifted vertically 
by equal amounts, (b) Angular average of the LDOS (top) 
and LDOS along nodal (♦) and anti-nodal (□) directions at 
the energies indicated by the arrows in (a). The strait lines 
show exponential behavior with a decay length A as indicated. 
Inset: LDOS anisotropy ratio (see text); smooth curves were 
obtained by interpolation of the lattice results. 

limit a(r) is maximum around r/a = 3, consistently with 
previous results The behavior of a(r) in the pseudo- 
gapped case is similar, although reduced by a factor ~ 25, 
making the maximum anisotropy smaller than 1%. 

Our model Cooperon propagator has deliberately few 
parameters. The good agreement between numerical re- 
sults and experimental data gives us confidence that it 
may contain the relevant input for interpreting vortex- 
core measurements in HTS. Our calculations show clearly 
that the same model and parameters used in the pseudo- 
gap regime also explain the properties of the core states 
in vortices. Whether the pseudogap part of the self- 
energy survives unchanged in the superconducting state 
or whether it reappears gradually only in the vortex cores 
should be considered an open problem, although our cal- 
culations in the homogeneous system and in the vortex, 
as well as some experimental evidence [|[ p8| , seem to 
favor the first term of the alternative. Moreover, the 
short-range part of the Cooperon propagator only repro- 
duces the experimental results if it is incoherent, i.e., 
not participating in the overall phase coherence below 
T c . This belongs today to the conventional wisdom, es- 
sentially because this part does not vanish in the center 
of vortices, and also because of the weak magnetic-field 
dependence of the pseudogap Finally, our model in- 
volves two different length scales, the superconducting 
coherence length and the range of the incoherent corre- 
lations. We note here that recent experiments by Pan et 
al. |2{J have shown the possible existence of two correla- 
tion lengths in BSCCO, although in their interpretation 
the shortest one is associated with the presence of the 
oxygen dopants in the BiO plane and is therefore not 
translationally invariant as in our approach. 
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According to our picture the pseudogap and super- 
conducting gap may have very different properties, and 
this seems to be supported by thermodynamic ]3(| and 
reflectivity |3l| experiments. The apparent contradiction 
with tunneling and photoemission data, which seem to 
show that the superconducting gap merges smoothly with 
the pseudogap at T c , can be resolved if the one-electron 
energy gap below T c turns out to be composed of both 
superconducting gap and pseudogap. This in fact hap- 
pens in our model, which interestingly produces a single 
gap structure in the DOS with a smooth evolution across 
T c , regardless of the relative values of A p and A S Q. 

Our theory being semi-phcnomcnological, it does not 
provide a deep understanding of the underlying micro- 
scopic phenomena. It is compatible, in principle, with 
any microscopic theory that would predict two energy 
scales — like the spin-charge separation scenario (which 
has the spinon gap and the superconducting gap as sep- 
arate entities) or the theory of stripe-induced supercon- 
ductivity (where, in one version, one-dimensional spinon 
field can exist above T c , the superconducting order be- 
ing created by Josephson coupling between the stripes) 
or the phase fluctuation picture, as emphasized below — 
and two length scales associated with these energy scales 
and appearing in the calculation of the Cooperon propa- 
gator. The picture that emerges for the pseudogap state 
in our model is analogous to a plane rotator model, where 
the two-dimensional spins would be formed at some tem- 
perature T* , but where the phase coherence only con- 
cerns the interaction between these spins (or preformed 
pairs), not the internal structure of the spins themselves. 
In the model, however, the two-dimensional spins appear 
only in the correlation function, not as entities located 
at definite sites. This picture, when proper temperature 
dependence is included in A S (T), could lead to a sec- 
ond cross-over temperature within the pseudogap state, 
perhaps seen in Nernst effect ]3^| , conductivity^], and 
Hall effect measurements [Q, temperature at which local 
phase coherence is established between neighboring spins, 
thus allowing the formation of fluctuating vortices. 

In conclusion, we display in this work a reasonably 
good agreement between a phenomenological theory of 
the pseudogap state and the experimental data for the 
vortex-core tunneling spectra in YBCO and BSCCO. Our 
results argue in favor of a common, pseudogap related, 
mechanism at the origin of the core states in these mate- 
rials. The incoherence and short-rangeness of the pseudo- 
gap correlations turn out to be the clue to the formation 
of these states, and endow them with properties (energy, 
amplitude, spatial decay) which are drastically different 
from the properties of conventional BCS core states. 

We wish to thank 0. Fischer, B. W. Hoogenboom, S. 
H. Pan, and J.-M. Triscone for valuable discussions. 
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